Computer evaluation of crop related industries

ABSTRACT

A process for optimizing a portfolio of products produced from a crop includes the use of an objective function to determine optimized quantities of the products in the portfolio. The objective function, for example, includes quantity terms for the products. The objective function may also include additional terms such as an energy term and/or a storage term. The energy term, for example, relates to an amount of energy required to produce the products. The storage term, for example, relates to the cost of storing products. The crop, for example, may be sugarcane, and the products, for example, may be sugar, molasses, bagasse, biofuel, electricity, and/or carbon credits.

RELATED APPLICATIONS

This application discloses subject matter similar to that disclosed in U.S. patent application Ser. No. 11/445,053 filed on May 31, 2006, U.S. patent application Ser. No. 11/872,826 filed on Oct. 16, 2007, and U.S. patent application Ser. No. 11/872,999 filed on Oct. 16, 2007.

TECHNICAL FIELD

The technical field of this disclosure relates, inter alia, to the evaluation of the trade offs relating to certain crop industries such as the sugar industry.

BACKGROUND

Many crop industries such as the sugar industry are becoming highly competitive due to globalization of competition. Cost effectiveness of a crop industry such as the sugar industry, for example, can be enhanced by optimally utilizing resources in the production of sugar products and by-products in different proportions in response to economic dynamics in the marketplace.

Traditionally, a crop industry such as the sugar industry has focused on the main crop product such as sugar. In the past, this strategy was justified as the other by-products (bio-fuels such as ethanol, butanol, bio-diesel or green-diesel, energy in form of electricity, bagasse, and molasses) were of secondary importance economically. In today's global economy, however, the price of the main crop product such as sugar fluctuates significantly as it is governed by global supply and demand. Also, because of the desire to decrease dependency on fossil fuels and because of an increased concern about global warming, the focus on and incentives for the production of bio-fuels have increased. The volatile prices of crude oil and electric power result in commensurate volatility in the price of bio-fuels. In addition, electric power generation using renewable sources (such as bagasse) fetches carbon credit revenues. For the sake of simplicity, the sugar industry will be considered as an example case of a crop industry.

The sugar industry has traditionally produced sugar from sugarcane juice. However, during the production of sugar, other major by-products such as molasses (the residual mass at the end of the crystallization stage) and bagasse (the fibrous residual after extraction of juices from sugarcane) are also produced. As indicated above, these by-products have traditionally been of less economical importance than sugar because of their low prices and demand.

Some sugar mills have converted molasses into bio-fuels by a fermentation process, and have converted bagasse into electricity using a co-generation method. However, because of low return on investment, the further conversion of by-products such as molasses and bagasse into sellable products like bio-fuels and electricity was not lucrative as many nations were mainly using crude oil (obtained by either domestic production or import) for their energy needs.

In today's global economy, where crude oil prices are threatening to become unaffordable, many countries are putting their efforts into developing energy self reliance as a part of their energy security strategies. Moreover, because of increased concern about global warming and stricter emission norms, both developed and developing countries are moving toward the adoption of environmentally friendly fuels. Hence, there is an increased demand for bio-fuels (fuel produced biologically) and green energy. Due to these changed circumstances, bio-fuel production and electricity generation offer additional sources of income to the sugar industry.

Sugar mills in today's changing world now have multiple economic routes for their resources. These mills can either divert the primary and/or secondary juice(s) to bio-fuel manufacturing by fermentation or they can use either or both juices for sugar production. Molasses can be either sold as animal feed and fertilizer or fermented to produce bio-fuels. Even bagasse can be burnt to produce electricity or sold as a raw material for the paper industry or fermented to produce bio-fuels. It should be noted that production of electricity is also possible from the spent wash generated during fermentation to produce bio-fuels such as ethanol. Moreover, the electricity generated using bagasse and spent wash is considered as green energy and earns carbon credits on excess energy (the energy above the internal needs of the sugar mill).

All of these economic alternatives require different time scales and processing costs and energies. In addition, government policies and regulations potentially affect the production levels of each of these end products.

Therefore, an agile decision support framework is needed. This framework should guide a crop industry such as the sugar industry to properly (i) evaluate these trade-offs under various what-if scenarios, and (ii) make optimal decisions related to the key intermediate resources, products, and by-products. Such a decision support system would help a crop industry such as a sugar mill to consistently drive its operations toward the optimization of profits and/or products and/or costs.

BRIEF DESCRIPTION OF THE DRAWINGS

Features, aspects, and advantages of the present invention will become better understood when the following detailed description is read with reference to the accompanying drawings in which:

FIG. 1 illustrates a process for optimizing user designed criterion in a crop industry such as a sugar mill;

FIG. 2 illustrates a computer that can be used to implement the process of FIG. 1; and,

FIGS. 3 and 4 provide additional detail for the process of FIG. 1; and,

FIG. 5 provides detailed information related to the application of the decision support system to the sugar industry.

DETAILED DESCRIPTION

An agile decision support system can help a crop industry to drive the industry in the most profitable manner. The decision support system described herein formulates and solves an optimization problem for a given what-if scenario. FIG. 4 illustrates steps involved in a decision support system 78. As shown in FIG. 4, at 80 in the decision support system 78, a user has the choice of selecting the best run criterion/criteria of the plant or mill. Thus, only the terms corresponding to the user selected criterion/criteria will be executed to provide the optimal solution. At 82, the decision support system 78 identifies the constraints, such as those described below, that influence the criterion/criteria selected at 80.

At 84, the decision support system 78 fetches and/or the user provides input data and parameters, such as those described below, related to the criterion/criteria selected at 80 and the constraints identified at 82.

At 86, the decision support system 78 determines the optimal outputs with respect to the criterion/criteria selected at 80. At 86, the decisions support system 78 can use so much of the objective function (47) as is pertinent to the selected constraints as well as to the input data and selected parameters. If the user changes the criterion/criteria as determined at 88, the new criterion/criteria are selected by the user at 90, and program flow returns to 82 to re-run the decision support program 78.

Accordingly, the user is assisted in simulating different what-if scenarios. The application of the decision support system 78 has been explained with reference sugar industry for the profit maximization scenario. In general, however, the decision support system 78 can be used to simulate any other what-if scenario in other crop industries such as where one of the product is bio-fuel and/or energy.

Sugar Industry Related Parameters

Sugarcane is the main raw material for a sugar mill. A typical sugar mill buys sugarcane from various farmers under contract with the sugar mill. Each sugar mill knows the planting date and variety of each crop of each different farmer. Accordingly, the sugar mill generates a harvest schedule for the entire harvest year. This harvest schedule though is revised over time. However, the quantity of sugarcane coming into the sugar mill on any given day is always between the maximum and minimum crushing capacity of the sugar mill.

The quality (generally represented by Pol in Cane or Recovery) of the sugarcane coming to the sugar mill on any given day is also known from the harvest schedule. Therefore, for technical simplicity, the amount of sugarcane harvested (QC_(d) ^(h)) on any given day d is equal to the amount of sugarcane crushed (QC_(d) ^(c)) on any given day d during the entire harvest season. After the sugarcane is crushed, the quantity of primary juice QPJ_(d) ^(p) obtained on day d and the quantity of secondary juice QSJ_(d) ^(p) obtained on day d are given by following equations: QPJ _(d) ^(p)=(PJC _(d))(QC _(d) ^(c))∀d  (1) QPJ _(d) ^(p)=(SJC _(d))(QC _(d) ^(c))∀d  (2) where PJC_(d) and SJC_(d) are the respective yields of the primary and secondary juices per ton of sugarcane crushed.

Sugarcane is generally crushed in two stages. The juice produced from the first stage is called primary juice, and the juice produced from the second stage is called secondary juice. Primary juice is more concentrated than secondary juice.

Though these yields depend on the quality of the cane crushed, the nominal values are, for example, thirty-seven tons of primary juice and sixty-five tons of secondary juice per 100 tons of cane crushed. Either these nominal fixed values for all harvest days or empirical relationships may be used to calculate the values for PJC_(d) and SJC_(d) for each given day, taking into consideration the quality of cane harvested on the given day. The empirical relationships are generally devised using domain knowledge, industry practice, and/or research findings, and are known to those skilled in the art. The sugar concentration CPJ_(d) in the primary juice for day d and the sugar concentration CSJ_(d) in the secondary juice for day d may be calculated using following set of equations:

$\begin{matrix} {{CPJ}_{d} = {\frac{\left( {POC}_{d} \right)\left( {PJSF}_{d} \right)}{{PJC}_{d}}\mspace{14mu}{\forall d}}} & (3) \\ {{CSJ}_{d} = {\frac{\left( {POC}_{d} \right)\left( {SJSF}_{d} \right)}{{SJC}_{d}}\mspace{14mu}{\forall d}}} & (4) \end{matrix}$ where POC_(d) is the quality Pol in the sugarcane on day d as predicted by the harvest schedule, PJSF_(d) is the fraction of the total sugar that is extracted from the sugarcane on day d for primary juice, and SJSF_(d) is the fraction of the total sugar that is extracted from the sugarcane on day d for secondary juice. The fractions PJSF_(d) and SJSF_(d) also depend on the quality of the sugarcane, and can be derived using domain knowledge, industry practice, and/or research findings. The nominal values for PJSF_(d) and SJSF_(d), for example, are in the ranges 0.7-0.8 and 0.3-0.2, respectively. The empirical relationships for these fractions are known to those skilled in the art.

Similarly, the quantity of bagasse produced on day d is known a priori taking into consideration the quantity and quality of sugarcane crushed on that day. Therefore, the quantity of produced bagasse is given by the following equation: QB _(d) ^(p)=(BC _(d))(QC _(d) ^(c))∀d  (5) where BC_(d) represents the yield of bagasse per ton of sugarcane. The nominal value for BC_(d) is, for example, 0.25 T. However, this value for BC_(d) depends on the quality of the sugarcane, and can be up to 0.4 T. The actual value for BC_(d) on day d can be calculated using the empirical relationship known to those skilled in art.

Sugar Industry Related Constraints

The sugar that is produced on day d depends on the split between primary and secondary juices that are used for sugar production. The fractions of the primary and secondary juices used for sugar production may be designated as fpj_(d) and fsj_(d), respectively. Hence, the quantity qs_(d) ^(p) of sugar produced on day d is given by the following equation: qs _(d) ^(p)=(CPJ _(d))(QPJ _(d) ^(p))(PJE _(d))(fpj _(d))+(CSJ _(d))(QSJ _(d) ^(P))(SJE _(d))(fsj _(d))∀d  (6) (∀d=for all d), where PJE_(d) and SJE_(d) are the mill efficiencies related to the conversion of the primary and secondary juices to sugar, respectively. Sugar is the main product from the sugar mill and is sold on the open market. However, the quantity of sugar to market qs_(d) ^(m) on day d cannot be more than the current stock of sugar. Hence, the following constraint applies:

$\begin{matrix} {{qs}_{d}^{m} \leq {{OPS} + {\sum\limits_{d^{\prime} = 1}^{d}{qs}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d - 1}{{qs}_{d^{\prime}}^{m}{\forall d}}}}} & (7) \end{matrix}$ where OPS represents the mill's opening stock of sugar at the start of the harvest season. The opening stock of sugar at the start of harvest season is the sugar that the sugar mill has on hand at the beginning of the harvest season.

The index d′ as used herein is an index notation for day count. If there is need to use more than one index notation for day count in same equation, then one of the index notations is designated with the prime.

Similarly on day d, the stock of sugar at the sugar mill should not be more than the capacity STCS_(d) of the sugar mill for storing sugar on day d. This constraint is given by the following equation:

$\begin{matrix} {{{OPS} + {\sum\limits_{d^{\prime} = 1}^{d}{qs}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}{qs}_{d^{\prime}}^{m}}} \leq {{STCS}_{d}\mspace{14mu}{\forall d}}} & (8) \end{matrix}$

It should be understood that producing sugar from primary juice is more economical compared to producing sugar from secondary juice. Hence, bio-fuels are produced from primary juice only after the entire supply of secondary juice is diverted for fermentation. Though primary and/or secondary juices can be used for the production of different bio-fuels such as ethanol, butanol, bio/green diesel, for discussion below will consider ethanol production as an example case. However, it is very straight forward for those skilled in the art to extend the stated formulation for any combination of bio-fuels.

The binary variable bpj_(d) may be used to indicate whether primary juice is to be used for ethanol production or not. Hence, if bpj_(d)=1, then primary juice is diverted to ethanol production, and if bpj_(d)=0, then primary juice is not diverted to ethanol production according to the following constraints: fsj _(d)≧(1−bpj _(d))∀d  (9) fpj _(d)≧(1−bpj _(d))∀d  (10)

It is assumed that the sugar mill has sufficient capacity to process entirely the primary and secondary juices. In other words, an fsj_(d)=1 and fpj_(d)=1 condition is allowed in the formulation without capacity constraints. However, when this assumption is not the case, then corresponding capacity constraints easily can be added.

Molasses Related Constraints

There are three types of molasses, namely A, B, and C molasses, that are produced during sugar production. The quantity qam_(d) ^(p) of type A molasses that is produced on day d, the quantity qbm_(d) ^(p) of type B molasses that is produced on day d, and the quantity qcm_(d) ^(p) of type C molasses that is produced on day d depend on what fractions of the primary and secondary juices are used for sugar production. These quantities are given by the following equations: qam _(d) ^(p)=(AMPJ _(d))(QPJ _(d) ^(p))(fpj _(d))+(AMSJ _(d))(QSJ _(d) ^(p))(fsj _(d))∀d  (11) qbm _(d) ^(p)=(BMPJ _(d))(QPJ _(d) ^(p))(fpj _(d))+(BMSJ _(d))(QSJ _(d) ^(p))(fsj _(d))∀_(d)  (12) qcm _(d) ^(p)=(CMPJ _(d))(QPJ _(d) ^(p))(fpj _(d))+(CMSJ _(d))(QSJ _(d) ^(p))(fsj _(d))∀d  (13) where AMPJ_(d) is the yield of type A molasses per ton of primary juice on day d, BMPJ_(d) is the yield of type B molasses per ton of primary juice on day d, CMPJ_(d) is the yield of type C molasses per ton of primary juice on day d, AMSJ_(d) is the yield of type A molasses per ton of secondary juice on day d, BMSJ_(d) is the yield of type B molasses per ton of secondary juice on day d, and CMSJ_(d) is the yield of type C molasses per ton of secondary juice on day d. These yields depend on the quality of the primary and secondary juices. The empirical relationship to calculate the values of these yields are known to those skilled in the art.

As discussed earlier, each type of molasses can be either sold to market or used for the production of bio-fuels. Though it is very much possible to produce different bio-fuels from molasses, for simplicity the discussion herein considers ethanol production as an example case. Modification of the formulation to consider any combination of bio-fuels is easy for those skilled in the art.

The total quantity allocated to market and used for ethanol manufacturing on day d should be less than the current storage as given by the following constraints:

$\begin{matrix} {{{qam}_{d}^{e} + {qam}_{d}^{m}} \leq {{OPAM} + {\sum\limits_{d^{\prime} = 1}^{d}{qam}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d - 1}{qam}_{d^{\prime}}^{e}} - {\sum\limits_{d^{\prime} = 1}^{d - 1}{{qam}_{d^{\prime}}^{m}\mspace{11mu}{\forall d}}}}} & (14) \\ {{{qbm}_{d}^{e} + {qbm}_{d}^{m}} \leq {{OPBM} + {\sum\limits_{d^{\prime} = 1}^{d}{qbm}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d - 1}{qbm}_{d^{\prime}}^{e}} - {\sum\limits_{d^{\prime} = 1}^{d - 1}{{qbm}_{d^{\prime}}^{m}\mspace{11mu}{\forall d}}}}} & (15) \\ {{{qcm}_{d}^{e} + {qcm}_{d}^{m}} \leq {{OPCM} + {\sum\limits_{d^{\prime} = 1}^{d}{qcm}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d - 1}{qcm}_{d^{\prime}}^{e}} - {\sum\limits_{d^{\prime} = 1}^{d - 1}{{qcm}_{d^{\prime}}^{m}\mspace{11mu}{\forall d}}}}} & (16) \end{matrix}$ where OPAM represents the opening stock for type A molasses on day d, OPBM represents the opening stock for type B molasses on day d, OPCM represents the opening stock for type C molasses on day d, qam_(d) ^(e) is the quantity of type A molasses that is used for ethanol production on day d, qbm_(d) ^(e) is the quantity of type B molasses that is used for ethanol production on day d, qcm_(d) ^(e) is the quantity of type C molasses that is used for ethanol production on day d, qam_(d) ^(m) is the quantity of type A molasses that is sold to market on day d, qbm_(d) ^(m) is the quantity of type B molasses that is sold to market on day d, and qcm_(d) ^(m) is the quantity of type C molasses that is sold to market on day d.

In addition, the quantity of each type of molasses available on day d within the sugar mill should be less than the storage capacity for each type as given by the following constraints:

$\begin{matrix} {{{OPAM} + {\sum\limits_{d^{\prime} = 1}^{d}{qam}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}{qam}_{d^{\prime}}^{e}} - {\sum\limits_{d^{\prime} = 1}^{d}{qam}_{d^{\prime}}^{m}}} \leq {{STCAM}_{d}\mspace{11mu}{\forall d}}} & (17) \\ {{{OPBM} + {\sum\limits_{d^{\prime} = 1}^{d}{qbm}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}{qbm}_{d^{\prime}}^{e}} - {\sum\limits_{d^{\prime} = 1}^{d}{qbm}_{d^{\prime}}^{m}}} \leq {{STCBM}_{d\;}{\forall d}}} & (18) \\ {{{OPCM} + {\sum\limits_{d^{\prime} = 1}^{d}{qcm}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}{qcm}_{d^{\prime}}^{e}} - {\sum\limits_{d^{\prime} = 1}^{d}{qcm}_{d^{\prime}}^{m}}} \leq {{STCCM}_{d\;}{\forall d}}} & (19) \end{matrix}$ where STCAM_(d) is the mill's storage capacity for molasses A on day d, STCBM_(d) is the mill's storage capacity for molasses B on day d, and STCCM_(d) is the mill's storage capacity for molasses C on day d.

Bagasse Related Constraints

The other by-product of sugar industry is bagasse which can also be sold to market or used for the production of bio-fuels. Though it is very much possible to produce different bio-fuels from bagasse, for simplicity the discussion herein considers ethanol production as an example case. Modification of the formulation to consider any combination of bio-fuels is easy for those skilled in the art. Additionally, bagasse can be burnt in co-generation facilities to produce electricity. However, its use and storage on day d is restricted by the available stock and storage capacity as given by the following constraints:

$\begin{matrix} {{{qb}_{d}^{e} + {qb}_{d}^{m} + {qb}_{d}^{el}} \leq {{OPB} + {\sum\limits_{d^{\prime} = 1}^{d}{QB}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d - 1}{\left( {{qb}_{d^{\prime}}^{e} + {qb}_{d^{\prime}}^{m} + {qb}_{d^{\prime}}^{el}} \right)\mspace{11mu}{\forall d}}}}} & (20) \\ {{{OPB} + {\sum\limits_{d^{\prime} = 1}^{d}{QB}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}\left( {{qb}_{d^{\prime}}^{e} + {qb}_{d^{\prime}}^{m} + {qb}_{d^{\prime}}^{el}} \right)}} \leq {{STCB}_{d}\mspace{11mu}{\forall d}}} & (21) \end{matrix}$ where qb_(d) ^(e) is the quantity of bagasse that is used for ethanol production on day d, qb_(d) ^(m) is the quantity of bagasse that is sold to market on day d, qb_(d) ^(el) is the quantity of bagasse that is burned for electricity on day d, OPB is the mill's opening stock of bagasse at the start of the harvest season, and STCB_(d) is the mill's storage capacity for bagasse on day d.

The quantity of bagasse allocated for burning in a co-generation facility should be less than the capacity of the co-generation facility as given by the following constraint: qb _(d) ^(el) ≦PCCP  (22) where PCCP is the processing capacity of a co-generation facility. Some co-generation plants need certain minimum quantities of bagasse to be burnt everyday to keep them operational without additional cost. The constraint related to this minimum quantity can be easily added into the formulation by those skilled in the art.

Bio-fuel Related Constraints

Bio-fuels are produced from primary juice, secondary juice, all types of molasses, and bagasse. However, for simplicity, this discussion consider ethanol production as an example case. Modification of the formulation described herein to consider any combination of bio-fuels is easy for those skilled in the art.

The yield of ethanol per ton from each of these above mentioned sources varies. Also, the time required to produce ethanol from these different sources is different. Thus, the quantity of ethanol that is produced is given by the following equation:

$\begin{matrix} {{qe}_{d}^{p} = {{({EAM})\left( {qam}_{d}^{e} \right)} + {({EBM})\left( {qbm}_{d}^{e} \right)} + {({ECM})\left( {qcm}_{d}^{e} \right)} + {({EB})\left( {qb}_{d}^{e} \right)} + {({EPJ})\left( {QPJ}_{d}^{p} \right)\left( {1 - {fpj}_{d}} \right)} + {({ESJ})\left( {QSJ}_{d}^{p} \right)\left( {1 - {fsj}_{d}} \right)\mspace{11mu}{\forall d}}}} & (23) \end{matrix}$ where qe_(d) ^(p) is the quantity of ethanol produced on day d, EAM is yield of ethanol per ton of type A molasses, EBM is yield of ethanol per ton of type B molasses, ECM is yield of ethanol per ton of type C molasses, EB is yield of ethanol per ton of bagasse, EPJ is yield of ethanol per ton of primary juice, and ESJ is yield of ethanol per ton of secondary juice,

The nominal values for ethanol yield from the different sources are known to those skilled in the art. For example, from one ton of type B molasses, 348 liters of ethanol is produced while the yield of ethanol per ton of type C molasses is around 250 liters. From one ton of raw secondary juice, approximately 48 liters of ethanol is produced. However, from one ton of heated secondary juice, the yield of ethanol is approximately 58 liters. In addition, the values of yield can depend on the quality of the source used. Hence, these ethanol yield values are dependent on the design of the ethanol plant and the quality of the source. Either the nominal value or the day dependent value, which can be obtained using the empirical relationships known to those skilled in the art, can be used.

The ethanol fermentation process can be either continuous or batch. The fermentation batch size and time in the case of batch operation or the source input rate to fermenter in the case of continuous operation depends on the type of source used. For example, the nominal values for time in the batch process, when secondary juice and types B and C molasses are used as the sources, are 15, 21 and 18 hours, respectively. In the continuous process, the total time required for all sources to feed to the ethanol fermenter(s) should be less than 1 day. Hence, the following constraint applies:

$\begin{matrix} {{\frac{{qam}_{d}^{e}}{FRAM} + \frac{{qbm}_{d}^{e}}{FRBM} + \frac{{qcm}_{d}^{e}}{FRCM} + \frac{{qb}_{d}^{e}}{FRB} + \frac{\left( {QPJ}_{d}^{p} \right)\left( {1 - {fpj}_{d}} \right)}{FRPJ} + \frac{\left( {QSJ}_{d}^{p} \right)\left( {1 - {fsj}_{d}} \right)}{FRSJ}} \leq {24\mspace{11mu}{\forall d}}} & (24) \end{matrix}$ where FRAM is the feed rate of type A molasses to an ethanol fermenter in a continuous process, FRBM is the feed rate of type B molasses to an ethanol fermenter in a continuous process, FRCM is the feed rate of type C molasses to an ethanol fermenter in a continuous process, FRB is the feed rate of bagasse to an ethanol fermenter in a continuous process, FRPJ is the feed rate of primary juice to an ethanol fermenter in a continuous process, and FRSJ is the feed rate of secondary juice to an ethanol fermenter in a continuous process.

In a batch operation, the quantities of the different sources allocated for ethanol production on day d can take more than one day time depending on the batch times. Hence, the following constraint applies:

$\begin{matrix} {{\sum\limits_{d^{\prime} = d}^{HD}\left( {\frac{\left( {qam}_{d^{\prime}}^{e} \right)({BTAM})}{BSAM} + \frac{\left( {qbm}_{d^{\prime}}^{e} \right)({BTBM})}{BSBM} + \frac{\left( {qcm}_{d^{\prime}}^{e} \right)({BTCM})}{BSCM} + \frac{\left( {qb}_{d^{\prime}}^{e} \right)({BTB})}{BSB} +} \right)}{\sum\limits_{d^{\prime} = d}^{HD}\left( {\frac{\left( {QPJ}_{d^{\prime}}^{p} \right)({BTPJ})\left( {1 - {fpj}_{d^{\prime}}} \right)}{BSPJ} + \frac{\left( {QSJ}_{d^{\prime}}^{p} \right)({BTSJ})\left( {1 - {fsj}_{d^{\prime}}} \right)}{BSSJ}} \right)} \leq {24\left( {{HD} - d} \right)\mspace{11mu}{\forall d}}} & (25) \end{matrix}$ where HD is the total number of harvest days, BETAM is the batch time for the ethanol fermentation of type A molasses, BTBM is the batch time for the ethanol fermentation of type B molasses, BTCM is the batch time for the ethanol fermentation of type C molasses, BTB is the batch time for the ethanol fermentation of bagasse, BTPJ is the batch time for the ethanol fermentation of primary juice, BTSJ is the batch time for the ethanol fermentation of secondary juice, BSPJ is the batch size of the primary juice used in ethanol fermentation, BSSJ is the batch size of the secondary juice used in ethanol fermentation, BSAM is the batch size of type A molasses used in ethanol fermentation, BSBM is the batch size of type B molasses used in ethanol fermentation, BSCM is the batch size of type C molasses used in ethanol fermentation, and BSB is the batch size of bagasse used in ethanol fermentation.

Ethanol has two grades, industry grade and fuel grade. The constraint given by equations (24) and (25) are for industry grade ethanol. Fuel grade ethanol is obtained from industry grade ethanol by a dehydration process. The processing capacity of the dehydration process may be designated as PCDP, and the quantity of industry grade ethanol sent to the dehydration process on day d to produce fuel grade ethanol may be designated as qe_(d) ^(fg). The following constraint then applies: qe _(d) ^(fg) ≦PCDP∀d  (26) The quantity qfge_(d) ^(p) of fuel grade ethanol that is produced on day d per liter of industry grade ethanol is approximately 0.96-0.97 liter and is given by the following equation: qfge _(d) ^(p)=(FGEIGE)(qe _(d) ^(fg))∀d  (27) where FGEIGE is the yield of fuel grade ethanol per ton of industry grade ethanol.

Both industry and fuel ethanol grades have market value. However, the quantity qe_(d) ^(m) of industry grade ethanol sold to the market on day d plus the quantity qe_(d) ^(fg) of industry grade ethanol that is used to produce fuel grade ethanol on day d should not be more than current available stock, as given by the following constraint:

$\begin{matrix} {{{qe}_{d}^{m} + {qe}_{d}^{fg}} \leq {{OPIGE} + {\sum\limits_{d^{\prime} = 1}^{d}{qe}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d - 1}{\left( {{qe}_{d^{\prime}}^{m} + {qe}_{d^{\prime}}^{fg}} \right)\mspace{11mu}{\forall d}}}}} & (28) \end{matrix}$ where OPIGE is the opening stock of industry grade ethanol at the start of the harvest season.

Similarly, the quantity qfge_(d) ^(m) of fuel grade ethanol sold to the market on day d should not be more than current available stock, as given by the following constraint:

$\begin{matrix} {{qfge}_{d}^{m} \leq {{OPFGE} + {\sum\limits_{d^{\prime} = 1}^{d}{qfge}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d - 1}{{qfge}_{d^{\prime}}^{m}\mspace{11mu}{\forall d}}}}} & (29) \end{matrix}$ where OPFGE is the opening stock of fuel grade ethanol at the start of the harvest season.

The quantity of industry grade ethanol available in the sugar mill on any given day d should not be more than the storage capacity STCIGE_(d) for industry grade ethanol on day d as given by the following constraint:

$\begin{matrix} {{{OPIGE} + {\sum\limits_{d^{\prime} = 1}^{d}{qe}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}\left( {{qe}_{d^{\prime}}^{m} + {qe}_{d^{\prime}}^{fg}} \right)}} \leq {{STCIGE}_{d}\mspace{11mu}{\forall d}}} & (30) \end{matrix}$ Similarly, the quantity of fuel grade ethanol available in the sugar mill on any given day d should not be more than the storage capacity STCFGE_(d) for fuel grade ethanol on day d as given by the following constraint:

$\begin{matrix} {{{OPFGE} + {\sum\limits_{d^{\prime} = 1}^{d}{qfge}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}{qfge}_{d^{\prime}}^{m}}} \leq {{STCFGE}_{d}\mspace{11mu}{\forall d}}} & (31) \end{matrix}$

During ethanol production, the sugar mill also produces spent wash as a byproduct. Approximately fourteen liters of spent wash is produced per liter of ethanol. However, this figure again depends on the type of source used and its quality. The total spent wash qsw_(d) ^(p) produced on any given day d is given by the following equation:

$\begin{matrix} {{qsw}_{d}^{p} = {{({SWAM})\left( {qam}_{d}^{e} \right)} + {({SWBM})\left( {qbm}_{d}^{e} \right)} + {({SWCM})\left( {qcm}_{d}^{e} \right)} + {({SWB})\left( {qb}_{d}^{e} \right)} + {({SWPJ})\left( {QPJ}_{d}^{p} \right)\left( {1 - {fpj}_{d}} \right)} + {({SWSJ})\left( {QSJ}_{d}^{p} \right)\left( {1 - {fsj}_{d}} \right){\forall d}}}} & (32) \end{matrix}$ where SWAM is the spent wash yield per ton of type A molasses, SWBM is the spent wash yield per ton of type B molasses, SWCM is the spent wash yield per ton of type C molasses, SWB is the spent wash yield per ton of bagasse, SWPJ is the spent wash yield per ton of primary juice, and SWSJ is the spent wash yield per ton of secondary juice.

The spent wash so generated is then converted into biogas and finally into electricity depending on the capacity availability.

Energy Related Constraints

All of the process operations such as juice evaporation, crystallization, ethanol fermentation and dehydration, etc. need the input of thermal, mechanical, and electrical energy. However, it is always possible for those skilled in the art to calculate the steam equivalent of the total energy requirement for all of these process operations. The total steam requirement sr_(d) for day d is given by the following equation:

$\begin{matrix} {{sr}_{d} = {{({SRCC})\left( {QC}_{d}^{c} \right)} + {\left( {SRPJS}_{d} \right)\left( {QPJ}_{d}^{p} \right)\left( {fpj}_{d} \right)} + {\left( {SRSJS}_{d} \right)\left( {QSJ}_{d}^{p} \right)\left( {fsj}_{d} \right)} + {\left( {SRPJE}_{d} \right)\left( {QPJ}_{d}^{p} \right)\left( {1 - {fpj}_{d}} \right)} + {\left( {SRSJE}_{d} \right)\left( {QSJ}_{d}^{p} \right)\left( {1 - {fsj}_{d}} \right)} + {({SRAME})\left( {qam}_{d}^{e} \right)} + {({SRBME})\left( {qbm}_{d}^{e} \right)} + {({SRCME})\left( {qcm}_{d}^{e} \right)} + {({SRBE})\left( {qb}_{d}^{e} \right)} + {({SRIGFGE})\left( {qe}_{d}^{fg} \right){\forall d}}}} & (33) \end{matrix}$ where SRCC is the steam requirement for crushing a ton of sugarcane, SRPJS_(d) is the steam requirement on day d to convert primary juice to sugar per ton of primary juice, SRSJS_(d) is the steam requirement on day d to convert secondary juice to sugar per ton of secondary juice, SRPJE_(d) is the steam requirement on day d to convert primary juice to ethanol per ton of primary juice, SRSJE_(d) is the steam requirement on day d to convert secondary juice to ethanol per ton of secondary juice, SRAME_(d) is the steam requirement on day d to convert type A molasses to ethanol per ton of type A molasses, SRBME_(d) is the steam requirement on day d to convert type B molasses to ethanol per ton of type B molasses, SRCME_(d) is the steam requirement on day d to convert type C molasses to ethanol per ton of type C molasses, SRBE_(d) is the steam requirement on day d to convert bagasse to ethanol per ton of bagasse, and SRIGFGE_(d) is the steam requirement on day d to convert industry grade ethanol to fuel grade ethanol per ton of industry grade ethanol.

In sugar mills, steam is generated using bagasse and/or fossil fuels according to the following equation: sp _(d)=(SGB)(qb _(d) ^(el))+(SGF)(qff _(d) ^(el))∀d  (34) where sp_(d) is the steam produced on day d, SGB is the steam economy from Bagasse, SGF is the steam economy from fossil fuel, and aff_(d) ^(el) is the quantity of fossil fuel to produce electricity on day d.

The total steam produced on any day d should be more than the requirement for steam on that given day. This constraint is important on days when sugarcane is not being crushed but ethanol fermentation and dehydration are running and is given by the following: sr _(d) ≦sp _(d)  (35)

Sugar mills can produce on day d more energy than they use and hence can send the surplus electricity to the power grid. This surplus of energy qel_(d) ^(p) is given by the following equation: qel _(d) ^(p)=(ELS)(sp _(d) −sr _(d))+(ELSW)(qsw _(d) ^(p))∀d  (36) where ELS is the electricity generated per ton of steam, and ELSW is the electricity generated per ton of spent wash.

It is assumed that the capacity to convert spent wash into electricity is sufficient. In addition, the steam generated can have a different pressure and temperature depending on the design of boiler used in the co-generation unit.

Demand Constraints

The correct way to obtain the optimal levels of various value added products is without any constraints on product quantities. However, sugar mills might have agreements with market players of value added products about the supplied quantities. The quantities subject to such agreements become demand constraints. It is safe to assume that these demand quantities are obtained taking into consideration processing capacity and supply agreements. Hence, the following set of constraints should be added to the decision support system formulation:

$\begin{matrix} \begin{matrix} {{\sum\limits_{d^{\prime} = 1}^{d}{qs}_{d^{\prime}}^{m}} \geq {\sum\limits_{d^{\prime} = 1}^{d}{DES}_{d^{\prime}}}} & {\forall d} \end{matrix} & (37) \\ \begin{matrix} {{\sum\limits_{d^{\prime} = 1}^{d}{qam}_{d^{\prime}}^{m}} \geq {\sum\limits_{d^{\prime} = 1}^{d}{DEAM}_{d^{\prime}}}} & {\forall d} \end{matrix} & (38) \\ \begin{matrix} {{\sum\limits_{d^{\prime} = 1}^{d}{qbm}_{d^{\prime}}^{m}} \geq {\sum\limits_{d^{\prime} = 1}^{d}{DEBM}_{d^{\prime}}}} & {\forall d} \end{matrix} & (39) \\ \begin{matrix} {{\sum\limits_{d^{\prime} = 1}^{d}{qcm}_{d^{\prime}}^{m}} \geq {\sum\limits_{d^{\prime} = 1}^{d}{DECM}_{d^{\prime}}}} & {\forall d} \end{matrix} & (40) \\ \begin{matrix} {{\sum\limits_{d^{\prime} = 1}^{d}{qb}_{d^{\prime}}^{m}} \geq {\sum\limits_{d^{\prime} = 1}^{d}{DEB}_{d^{\prime}}}} & {\forall d} \end{matrix} & (41) \\ \begin{matrix} {{\sum\limits_{d^{\prime} = 1}^{d}{qe}_{d^{\prime}}^{m}} \geq {\sum\limits_{d^{\prime} = 1}^{d}{DEIGE}_{d^{\prime}}}} & {\forall d} \end{matrix} & (42) \\ \begin{matrix} {{\sum\limits_{d^{\prime} = 1}^{d}{qfge}_{d^{\prime}}^{m}} \geq {\sum\limits_{d^{\prime} = 1}^{d}{DEFGE}_{d^{\prime}}}} & {\forall d} \end{matrix} & (43) \\ \begin{matrix} {{qel}_{d}^{p} \geq {DEEL}_{d}} & {\forall d} \end{matrix} & (44) \end{matrix}$ where DES_(d) is the demand for sugar on day d, DEAM_(d) is the demand for type A molasses on day d, DEBM_(d) is the demand for type B molasses on day d, DECM_(d) is the demand for type C molasses on day d, DEB_(d) is the demand for bagasse on day d, DEIGE_(d) is the demand for industry grade ethanol on day d, DEFGE_(d) is the demand for fuel grade ethanol on day d, and DEEL_(d) is the demand for electricity on day d.

The production of any value added product till day d can be more than the given demand on any given day d.

Carbon Credits Related Constraints

The electricity generated using bagasse is considered green energy and earns carbon credits on the excess (surplus) energy that is generated as indicated by the following equation:

$\begin{matrix} \begin{matrix} {{ccbel}_{d} = {({ELS})\left( {1 + \frac{TL}{100}} \right)\left( {{NCECR} - {NCEBR}} \right)\left( {{sp}_{d} - {sr}_{d}} \right)}} & {\forall d} \end{matrix} & (45) \end{matrix}$ where ccbel_(d) is the carbon credits that are earned on day d from electricity that is generated using bagasse, TL is the electricity transmission loss, NCECR is the net carbon emission for electricity generated using coal, and NCEBR is the net carbon emission for electricity generated using bagasse.

The electricity generated using spent wash is also considered green energy and earns carbon credits on the excess (surplus) energy that is generated as indicated by the following equation:

$\begin{matrix} \begin{matrix} {{ccswel}_{d} = {({ELSW})\left( {1 + \frac{TL}{100}} \right)\left( {{NCECR} - {NCESWR}} \right)\left( {qsw}_{d}^{p} \right)}} & {\forall d} \end{matrix} & (46) \end{matrix}$ where ccswel_(d) is the carbon credits that are earned on day d from electricity that is generated using spent wash, and NCESWR is the net carbon emission for electricity generated using spent wash.

The transmission losses given by TL vary from country to country and vicinity of the area using the electricity generated. These transmission losses generally vary in the range of 20-40%. Nominal values for the net carbon emission using coal is, for example, 950 g/KWH, and using bagasse is, for example, 120 g/KWH (bagasse route). The nominal value for the net carbon emission using spent wash is generally less than that using bagasse. The values of NCEBR and NCESWR depend on farming practices and use of locomotives, tractors, etc. during farming and can change from location to location. The bio-fuels generated using sugar industry by-products is also a form of green energy. In today's carbon credit system, bio-fuels do not earn any credit. However, if bio-fuels earn credit in the future, constraints to calculate the carbon credit from bio-fuels can be easily added.

Objective Function

One objective herein is to determine the optimal production levels for value added products taking into consideration the various economical trade-offs. The objective function is then used to maximize profit over the entire harvest season.

An example of an objective function to determine these optimum production levels is given by the following equation:

$\begin{matrix} {{\max\mspace{11mu}{OBJ}} = {{\sum\limits_{d}\left( {{\left( {PS}_{d} \right)\left( {qs}_{d}^{m} \right)} + {\left( {PIGE}_{d} \right)\left( {qe}_{d}^{m} \right)} + {\left( {PFGE}_{d} \right)\left( {qfge}_{d}^{m} \right)} + {\left( {PB}_{d} \right)\left( {qb}_{d}^{m} \right)}} \right)} + {\sum\limits_{d}\left( {{\left( {PAM}_{d} \right)\left( {qam}_{d}^{m} \right)} + {\left( {PBM}_{d} \right)\left( {qbm}_{d}^{m} \right)} + {\left( {PCM}_{d} \right)\left( {qcm}_{d}^{m} \right)}} \right)} + {\sum\limits_{d}\left( {{\left( {PEU}_{d} \right)\left( {qel}_{d}^{m} \right)} + {\left( {PCC}_{d} \right)\left( {{ccbel}_{d} - {ccswel}_{d}} \right)} + {\left( {PFF}_{d} \right)\left( {qff}_{d}^{el} \right)}} \right)} - {\sum\limits_{d}\left\lbrack {{STPSF}_{d} + {\left( {STPSV}_{d} \right)\left( {{OPS} + {\sum\limits_{d^{\prime} = 1}^{d}{qs}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}{qs}_{d^{\prime}}^{m}}} \right)}} \right\rbrack} - {\sum\limits_{d}\left\lbrack {{STPIGEF}_{d} + {\left( {STPIGEV}_{d} \right)\left( {{OPIGE} + {\sum\limits_{d^{\prime} = 1}^{d}{qe}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}\left( {{qe}_{d^{\prime}}^{m} + {qe}_{d^{\prime}}^{fg}} \right)}} \right)}} \right\rbrack} - {\sum\limits_{d}\left\lbrack {{STPFGEF}_{d} + {\left( {STPFGEV}_{d} \right)\left( {{OPFGE} + {\sum\limits_{d^{\prime} = 1}^{d}{qfge}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}{qfge}_{d^{\prime}}^{m}}} \right)}} \right\rbrack} - {\sum\limits_{d}\left\lbrack {{STPAMF}_{d} + {\left( {STPAMV}_{d} \right)\left( {{OPAM} + {\sum\limits_{d^{\prime} = 1}^{d}{qam}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}{qam}_{d^{\prime}}^{e}} - {\sum\limits_{d^{\prime} = 1}^{d}{qam}_{d^{\prime}}^{m}}} \right)}} \right\rbrack} - {\sum\limits_{d}\left\lbrack {{STPBMF}_{d} + {\left( {STPBMV}_{d} \right)\left( {{OPBM} + {\sum\limits_{d^{\prime} = 1}^{d}{qbm}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}{qbm}_{d^{\prime}}^{e}} - {\sum\limits_{d^{\prime} = 1}^{d}{qbm}_{d^{\prime}}^{m}}} \right)}} \right\rbrack} - {\sum\limits_{d}\left\lbrack {{STPCMF}_{d} + {\left( {STPCMV}_{d} \right)\left( {{OPCM} + {\sum\limits_{d^{\prime} = 1}^{d}{qcm}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}{qcm}_{d^{\prime}}^{e}} - {\sum\limits_{d^{\prime} = 1}^{d}{qcm}_{d^{\prime}}^{m}}} \right)}} \right\rbrack} - {\sum\limits_{d}{\left\lbrack {{STPBF}_{d} + {\left( {STPBV}_{d} \right)\left( {{OPB} + {\sum\limits_{d^{\prime} = 1}^{d}{QB}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}\left( {{qb}_{d^{\prime}}^{m} + {qb}_{d^{\prime}}^{el} + {qb}_{d^{\prime}}^{e}} \right)}} \right)}} \right\rbrack{\forall d}}}}} & (47) \end{matrix}$ where PS_(d) is the price of sugar on day d, PIGE_(d) is the price of industry grade ethanol on day d, PFGE_(d) is the price of fuel grade ethanol on day d, PB_(d) is the price of bagasse on day d, PAM_(d) is the price of type A molasses on day d, PBM_(d) is the price of type B molasses on day d, PCM_(d) is the price of type C molasses on day d, PEU_(d) is the price of electricity on day d, PCC_(d) is the price of carbon credits on day d, PFF_(d) is the price of fossil fuel on day d, STPSF_(d) is the fixed price for storage of sugar on day d, STPSV_(d) is the variable price for storage of sugar on day d, STPIGEF_(d) is the fixed price for storage of industry grade ethanol on day d, STPIGEV_(d) is the variable price for storage of industry grade ethanol on day d, STPFGEF_(d) is the fixed price for storage of fuel grade ethanol on day d, STPFGEV_(d) is the variable price for storage of fuel grade ethanol on day d, STPAMF_(d) is the fixed price for storage of type A molasses on day d, STPAMV_(d) is the variable price for storage of type A molasses on day d, STPBMF_(d) is the fixed price for storage of type B molasses on day d, STPBMV_(d) is the variable price for storage of type B molasses on day d, STPCMF_(d) is the fixed price for storage of type C molasses on day d, STPCMV_(d) is the variable price for storage of type C molasses on day d, STPBF_(d) is the fixed price for storage of bagasse on day d, and STPBV_(d) is the variable price for storage of bagasse on day d.

Quantities represented with lower case letters herein are variables, and quantities represented with upper case lettering herein are parameters. The lower case quantities in equation (47) are the variables of the objective function. The values of these variables that produce a maximum output of equation (47) are to be determined from the objective function. The values of the parameters in equation (47) are inputs or constants.

When the objective function of equation (47) is used to determine the best mix of sugar products and byproducts at a future day d, the parameters in equation (47) will need to be forecasted. Therefore, the accuracy in determining the mix of sugar products and byproducts that produce maximum profits will depend on the accuracy of the forecasts.

The first three lines of equation (47) represent the revenues that can be derived from various sugar products and byproducts. Subtracted from these revenues are the storage costs for these products and byproducts represented by the remaining seven lines of equation (47).

It should be understood that the objective function of equation (47) can be simplified or modified depending upon which value added products are to be produced and which constraints are to be considered.

The formulation given by equation (47) is a mixed integer linear programming formulation. The stated objective addresses the various trade-offs to determine optimal levels of production for different value added sugar related products and byproducts. Also, the formulation can be used to analyze the incremental benefits due to an increase in fermentation, co-generation or storage capacities.

Because the formulation given by equation (47) is a mixed integer linear programming formulation, an iterative mixed integer linear program can be applied to process the objective function. Instead of using iterative mixed integer linear programming, non-iterative mixed integer non-linear programming may be used.

FIG. 1 is a flow chart of a process 10 for determining optimized quantities for (i) the amount of sugar to be sold to market on day d (qs_(d) ^(m)), (ii) the amount of industry grade ethanol to be sold to market on day d (qe_(d) ^(m)), (iii) the amount of fuel grade ethanol to be sold to market on day d (qfge_(d) ^(m)), (iv) the amount of bagasse to be sold to market on day d (qb_(d) ^(m)), (v) the amount of type A molasses to be sold to market on day d (qam_(d) ^(m)), (vi) the amount of type B molasses to be sold to market on day d (qbm_(d) ^(m)), (vii) the amount of type C molasses to be sold to market on day d (qcm_(d) ^(m)), (viii) the amount of electricity to be produced on day d (qel_(d) ^(p)), (ix) the number of carbon credits to be received from electricity generated on day d by use of bagasse (ccbel_(d)), (x) the number of carbon credits to be received from electricity generated on day d by use of spent wash (ccswel_(d)), and (xi) the amount of fossil fuel used to generate electricity on day d (aff_(d) ^(el)).

The process 10 also determines optimized stored quantities for (i) the amount of sugar that is produced but not sold to market

$\begin{matrix} {\left( {{\sum\limits_{d^{\prime} = 1}^{d}{qs}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}{qs}_{d^{\prime}}^{m}}} \right),} & ({ii}) \end{matrix}$ the amount of industry grade ethanol that is produced but not sold to market nor used for the production of fuel grade ethanol

$\begin{matrix} {\left( {{\sum\limits_{d^{\prime} = 1}^{d}{qe}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}\left( {{qe}_{d^{\prime}}^{m} + {qe}_{d^{\prime}}^{fg}} \right)}} \right),} & ({iii}) \end{matrix}$ the amount of fuel grade ethanol that is produced but not sold to market

$\begin{matrix} {\left( {{\sum\limits_{d^{\prime} = 1}^{d}{qfge}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}{qfge}_{d^{\prime}}^{m}}} \right),} & ({iv}) \end{matrix}$ the amount of type A molasses that is produced but not sold to market nor used in ethanol production

$\begin{matrix} {\left( {{\sum\limits_{d^{\prime} = 1}^{d}{qam}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}{qam}_{d^{\prime}}^{e}} - {\sum\limits_{d^{\prime} = 1}^{d}{qam}_{d^{\prime}}^{m}}} \right),} & (v) \end{matrix}$ the amount of type B molasses that is produced but not sold to market nor used in ethanol production

$\begin{matrix} {\left( {{\sum\limits_{d^{\prime} = 1}^{d}{qbm}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}{qbm}_{d^{\prime}}^{e}} - {\sum\limits_{d^{\prime} = 1}^{d}{qbm}_{d^{\prime}}^{m}}} \right),} & ({vi}) \end{matrix}$ the amount of type C molasses that is produced but not sold to market nor used in ethanol production

$\begin{matrix} {\left( {{\sum\limits_{d^{\prime} = 1}^{d}{qcm}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}{qcm}_{d^{\prime}}^{e}} - {\sum\limits_{d^{\prime} = 1}^{d}{qcm}_{d^{\prime}}^{m}}} \right),{and}} & ({vii}) \end{matrix}$ the amount of bagasse that is produced but not sold to market nor used in ethanol production nor used for the generation of electricity

$\left( {{\sum\limits_{d^{\prime} = 1}^{d}{QB}_{d^{\prime}}^{p}} - {\sum\limits_{d^{\prime} = 1}^{d}\left( {{qb}_{d^{\prime}}^{m} + {qb}_{d^{\prime}}^{el} + {qb}_{d^{\prime}}^{e}} \right)}} \right).$

Accordingly, as shown in FIG. 1, the relevant input data, which is used in the equations and constraints (1)-(47) and which is maintained in a database 12, is loaded from memory at 14, and the objective function of equation (47) is maximized at 16 using the equations and constraints (1)-(46) as provided at 18.

The optimized quantities that maximize the objective function are provided at 20 as an output of the process 10. This output can be a display, such as on a monitor or in printed form. Alternatively or additionally, this output could be orders and/or contracts to sell and/or to produce the optimized quantities. Other forms of output could also be provided. Thus, the output 20 represents the decisions that can be made as to how the sugar mill can be optimally run.

The process corresponding to the flow chart of FIG. 1 can be executed in connection with a computer 30 shown in FIG. 2. The computer 30 includes a processor 32, a memory 34, an input device(s) 36, and an output device(s) 38.

The input device(s) 36 may be a mouse, a keyboard, etc. capable of inputting data to the processor 32. The input device(s) 36 may be used to input the data described above. All of this data may be stored in the memory 34.

The output device(s) 38 may be a monitor, a printer, etc. capable of outputting the quantities discussed above.

The memory 34 stores the input data, the modeling procedure 10 shown in FIG. 2, and/or the outputs provided by the output device(s) 38.

FIG. 3 explains the information flow from a plantation and harvest scheduler 42 to the decision support system 10. The information obtained from the plantation and harvest scheduler 42, includes for example, the amount of sugarcane harvested (QC_(d) ^(h)) on any given day d, the amount of sugarcane crushed (QC_(d) ^(c)) on any given day d during the entire harvest season, the respective yields of primary and secondary juices PJC_(d) and SJC_(d), the sugar concentration CPJ_(d) in the primary juice for day d, the sugar concentration CSJ_(d) in secondary juice for day d, the quality in the sugarcane on day d as predicted by the harvest schedule (POC_(d)), the fraction of the total sugar that is extracted from the sugarcane on day d for primary juice (PJSF_(d)), and the fraction of the total sugar that is extracted from the sugarcane on day d for secondary juice (SJSF_(d)), and the yield of bagasse per ton of sugarcane (BC_(d)). The parameters discussed herein can be provided by sources other than the plantation and harvest scheduler 42.

The parameters required for decision support system optimization are calculated by a parameter calculation module 44. For example, the parameter calculation module 44 contains the empirical relationships to calculate the values of the different parameters used herein. These empirical relationships are known to those skilled in the art.

FIG. 4 illustrates another version of the program of FIG. 1. This program version is designated 78 in FIG. 4. As shown in FIG. 4, a user has the choice of selecting the best run criterion/criteria of the plant or mill. Thus, only that part of equation (47) corresponding to the user selected criterion/criteria will be executed to provide the optimal solution. For example, that user may select to optimize the quantity of sugar to be sold to market on day d (qs_(d) ^(m)), the quantity of bagasse to be sold to market on day d (qb_(d) ^(m)), and the quantity of type A molasses to be sold to market on day d (qam_(d) ^(m)). The user makes this selection at 80 of the program shown in FIG. 4. Accordingly, only the following part of equation (47) is used:

${\max\mspace{11mu}{OBJ}} = {{\sum\limits_{d}\left( {{\left( {PS}_{d} \right)\left( {qs}_{d}^{m} \right)} + {\left( {PB}_{d} \right)\left( {qb}_{d}^{m} \right)}} \right)} + {\sum\limits_{d}{\left( {PAM}_{d} \right)\left( {qam}_{d}^{m} \right)}}}$

At 82, the program 78 will determine the constraints that influence the criterion/criteria selected at 80. In the example above, the program at 82 determines the constraints as given by Equations 1-6, 7, 8, 11, 14, 20, 21, 37, 38, and 41 as well as the following constraints: fpj _(d)=1∀d fsj _(d)=1∀d qam _(d) ^(e)=0∀d qb _(d) ^(e)=0∀d qb _(d) ^(el)=0∀d

The user will supply input data and parameters related to the criterion/criteria selected at 80, and at 84 the program 78 fetches these inputs and parameters related to the constraints determined at 82.

At 86, the program 78 determines the optimal outputs with respect to the criterion/criteria selected at 80. Further to the above example, the program 78 at 80 determines the optimized quantity of sugar to be sold to market on day d (qs_(d) ^(m)), the optimized quantity of bagasse to be sold to market on day d (qb_(d) ^(m)), and the optimized quantity of type A molasses to be sold to market on day d (qam_(d) ^(m)). These optimized quantities are also displayed to the user at 86.

If the user changes the criterion/criteria as determined at 88, the new criterion/criteria are selected by the user at 90, and program flow returns to 82 to re-run the program 78.

Accordingly, the user is assisted in simulating different what-if scenarios.

FIG. 5 shows the input information that is required by the process 10 with the sugar industry as an example case. This input information includes demand data 50. The demand data 50 includes the sugar demand profile DES_(d), the molasses demand profiles DEAM_(d), DEBM_(d), and DECM_(d), the bagasse demand profile DEB_(d), the bio-fuel demand profiles such as DEIGE_(d) and DEFGE_(d), and the power demand profile DEEL_(d).

This input information also includes pricing data 52. The pricing data 52 includes the forecasted sugar price, the forecasted bio-fuel price, the forecasted bagasse price, the forecasted molasses price, the forecasted power price, the forecasted molasses storage cost, the forecasted bagasse storage cost, the forecasted sugar storage cost, and the forecasted bio-fuel storage cost.

The input information additionally includes capacity data 54. The capacity data 54 includes fermentation plant capacity PCDP, cogeneration plant capacities PCCP, molasses storage capacities STCAM_(d), STCBM_(d), and STCCM_(d), bagasse storage capacity STCB_(d), sugar storage capacity STCS_(d), and bio-fuel storage capacities STCIGE_(d) and STCFGE_(d).

This input information also includes raw material data 56. The raw material data 56 includes the quantity of sugarcane to be crushed (QC_(d) ^(c)), and the POC (pure obtainable sugar) in the sugarcane to be crushed.

The input information further includes limitations 58 imposed by the policies of any relevant governments. The limitations 58 includes sugar production limitations, bio-fuel production limitations, export limitations, and sugar inventory limitations. The constraints related to government mandates can be easily added by those skilled in the art.

The input information also includes parameter values 60. The parameter values 60 include primary and secondary juice yields PJC_(d) and PSC_(d), bagasse yield BC_(d), molasses yields from primary and secondary juices AMPJ_(d), AMSJ_(d), BMPJ_(d), BMSJ_(d), CMPJ_(d), and CMSJ_(d), bio-fuel yields from juice, molasses, and bagasse such as EPJ, ESJ, EAM, EBM, ECM, EB, and FGEIGE, percentages of POC in primary and secondary juices represented by PJSF_(d) and SJSF_(d), and feed rates for juice, molasses, and bagasse for bio-fuel production such as FRPJ, FRSJ, FRAM, FRBM, FRCM, and FRB.

The input data include any other parameters as described herein or otherwise.

The block 20 of FIG. 1 is also represented by the blocks 62, 64, 66, 68, and 70 of FIG. 5. As shown in FIG. 5, the block 62 provides the optimized sugar profiles which include production and sale of sugar on day d, the block 64 provides the optimized bio-fuel profiles which include production and sale of industry and fuel grade bio-fuel on day d, the block 66 provides the optimized molasses (A, B, and C) profiles which include production, utilization, and sale of different (A, B and C) molasses on day d, the block 68 provides the optimized power that is sold to market on day d, and the block 70 provides the optimized bagasse profile which includes production, utilization, and sale of bagasse on day d.

Certain modifications of the present invention have been discussed above. Other modifications of the present invention will occur to those practicing in the art of the present invention. For example, the present invention has been described above in connection with sugarcane, sugar beet, and sugar related products. However, the present invention could be used in connection with other crops where one of the product is bio-fuel and/or energy. Typical examples of such other crops include wheat, sweet sorghum, corn, etc.

As another example, the objective function of equation (47) might rely on fewer products and therefore fewer storage requirements for those fewer products.

The framework described herein can be used to obtain decisions as to how to best run a mill by taking into consideration known data/information. There are three aspects to this framework: 1) the known data or information; 2) the best run of the mill (based on some criteria); and, 3) decisions such as the output 20.

The known data/information includes, but is not limited to, a) the constraints related to mill capacities, b) pricing information of the products and raw materials, c) government mandates related to products and emissions, d) demand targets of products, and e) the energy requirement for product production.

The best run of the mill is based on maximization or minimization of any combination of the following criteria but is not limited to: a) product quantities to market; b) product revenue; c) carbon credits; d) storage cost of products; e) storage quantities of products; and, f) energy generated using products

The decisions (outputs) to achieve the best run can be, but are not limited to, a) quantities of products produced and/or sold in market; b) storage quantities of products; c) quantity of energy and/or electricity produced using products; and, d) quantity of electricity to be imported/exported.

Accordingly, the description of the present invention is to be construed as illustrative only and is for the purpose of teaching those skilled in the art the best mode of carrying out the invention. The details may be varied substantially without departing from the spirit of the invention, and the exclusive use of all modifications which are within the scope of the appended claims is reserved. 

What is claimed is:
 1. A computer implemented method for optimizing production of products in a product portfolio of a crop industry comprising: computer receiving user selected variables for products to be produced for the crop industry product portfolio to be optimized; computer determining constraints relating to the products to be produced for the crop industry product portfolio to be optimized and to the variables selected by the user; computer fetching input parameters related to the determined constraints and to an objective function to be optimized; and, computer optimizing quantities of the products to be produced and supplied to market based on the constraints and the objective function, wherein the objective function includes terms related to the parameters and the user selected variables.
 2. The method of claim 1 wherein the optimizing of quantities comprises maximizing the objective function.
 3. The method of claim 1 wherein the variables include quantity terms and an energy term, wherein the quantity terms relate to the quantities of the products to be optimized, and wherein the energy term includes an amount of energy required to produce the products.
 4. The method of claim 3 wherein the energy in the energy term also includes an amount of energy produced from the products.
 5. The method of claim 1 wherein the variables include quantity terms and storage terms, wherein the quantity terms relate to the quantities of the products to be optimized, and wherein the storage terms includes storage costs for storing the products that are produced but not supplied to market.
 6. The method of claim 1 wherein the variables include quantity terms and a carbon credit term, wherein the quantity terms relate to the quantities of the products to be optimized, and wherein the carbon credit term relates to carbon credits realizable from production of electricity from the products.
 7. The method of claim 6 wherein the carbon credit term further relates to fossil fuels required to produce the products.
 8. The method of claim 1 wherein the variables include quantity terms, an energy term, and storage terms, wherein the quantity terms relate to the quantities of the products to be optimized, wherein the energy term includes an amount of energy required to produce the products, and wherein the storage terms includes storage costs for storing the products that are produced but not supplied to market.
 9. The method of claim 8 wherein the energy in the energy term also includes an amount of energy produced from the products.
 10. The method of claim 1 wherein the variables include quantity terms, an energy term, and a carbon credit term, wherein the quantity terms relate to the quantities of the products to be optimized, wherein the energy term includes an amount of energy required to produce the products, and wherein the carbon credit term relates to carbon credits realizable from production of electricity from the products.
 11. The method of claim 10 wherein the energy in the energy term also includes an amount of energy produced from the products, and wherein the carbon credit term further relates to fossil fuels required to product the products.
 12. The method of claim 1 wherein the variables include quantity terms, storage terms, and a carbon credit term, wherein the quantity terms relate to the quantities of the products to be optimized, wherein the storage terms include storage costs for storing the products that are produced but not supplied to market, and wherein the carbon credit term relates to carbon credits realizable from production of electricity from the products.
 13. The method of claim 12 wherein the carbon credit term further relates to fossil fuels required to product the products.
 14. The method of claim 1 wherein the variables include quantity terms, an energy term, storage terms, and a carbon credit term, wherein the quantity terms relate to the quantities of the products to be optimized, wherein the energy term includes an amount of energy required to produce the products, wherein the storage terms include storage costs for storing the products that are produced but not supplied to market, and wherein the carbon credit term relates to carbon credits realizable from production of electricity from the products.
 15. The method of claim 14 wherein the energy in the energy term also includes an amount of energy produced from the products, and wherein the carbon credit term further relates to fossil fuels required to product the products.
 16. The method of claim 1 wherein the products include sugar.
 17. The method of claim 16 wherein the products further include electricity.
 18. The method of claim 16 wherein the products further include bio-fuel.
 19. The method of claim 18 wherein the bio-fuel includes industry grade bio-fuel and fuel grade bio-fuel.
 20. The method of claim 16 wherein the products further include carbon credits.
 21. The method of claim 16 wherein the products further include molasses.
 22. The method of claim 21 wherein the molasses include type A molasses, type B molasses, and type C molasses.
 23. The method of claim 16 wherein the products further include bagasse.
 24. The method of claim 1 wherein the variables include quantity terms, wherein the quantity terms are based on the quantities of the products to be optimized, and wherein the quantity terms include corresponding market supply prices for the products.
 25. The method of claim 1 wherein the optimizing of quantities includes providing decisions as to the products.
 26. The method of claim 25 wherein the decisions include decisions on quantities of the products produced, sold, and/or stored, quantities of energy produced, and/or quantities of energy to imported and/or exported.
 27. The method of claim 1 wherein the constraints include pricing information, government mandates, demand profiles, plant capacity limitation, and/or energy requirements.
 28. The method of claim 1 wherein the optimized quantities include product revenue, carbon credits, storage costs, stored quantities, and/or generated energy.
 29. The method of claim 1 wherein the optimizing of quantities comprises iteratively solving the objective function.
 30. The method of claim 1 wherein the crop industry is such that one of the products includes a bio-fuel or a bio-energy.
 31. The method of claim 30 wherein the biofuel includes bio-ethanol, bio-butanol, bio-hydrogen, bio-gas, bio-diesel, green diesel, and/or any combination thereof.
 32. The method of claim 30 wherein the bio-energy includes thermal energy, mechanical energy, electrical energy, and/or any combination thereof. 